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ABSTRACT 

We consider the interaction of differential rotation and toroidal fields that are current- 
free in the gap between two corotating axially unbounded cylinders. It is shown that 
nonaxisymmetric perturbations are unstable if the rotation rate and Alfven frequency 
of the field are of the same order almost independent of the magnetic Prandtl number 
Pm. For the very steep rotation law fl cx R~ 2 (the Rayleigh limit) this Azimuthal 
MagnetoRotational Instability (AMRI) scales with the ordinary Reynolds number and 
the Hartmann number, which allows a laboratory experiment with liquid metals like 
sodium or gallium in a Taylor- Couette container. The growth rate of AMRI scales with 
fl 2 in the low-conductivity limit and with f2 in the high-conductivity limit. For the 
weakly nonlinear system the numerical values of the kinetic energy and the magnetic 
energy are derived for magnetic Prandtl numbers between 0.05 and unity. We find that 
the magnetic energy scales with the magnetic Reynolds number Rm, while the kinetic 
energy scales with Rm/ VPm. The resulting turbulent Schmidt number as the ratio of 
the diffusion coefficients of angular momentum and a passive scalar (such as lithium) 
is of order 20 for Pm = 1, but for small Pm drops to order unity. Hence, in a stellar 
core with fossil fields and steep rotation law the transport of angular momentum by 
AMRI is always accompanied by an intense mixing of the plasma, until the rotation 
becomes rigid. 

Key words: magnetohydrodynamics (MHD) - instabilities - stars: magnetic fields 
- stars: interiors. 



1 INTRODUCTION 

A hydrodynamically stable rotation law may become unsta- 
ble under the presence of a sufficiently strong uniform axial 
field. For this so-called magnetorotational instability (MRI) 
the axial field gives the catalyst which makes the rotation 
law unstable without any modification of the uniform field 
(Velikhov 1959). Of course, by this mechanism a protoplan- 
etary disk with a Keplerian rotation law should also become 
unstable if it is not too cold - or (almost) equivalently - 
if the field is strong enough. For an estimation of the criti- 
cal magnetic field strength it is enough to use the condition 
S > 1, where the Lundquist number 



B,H 



(1) 



with B z the amplitude of the axial field, H the semi- 
thickness of the disk, and r\ the microscopic magnetic difhi- 
sivity. Plotted in the S-Rm plane, with 



Rm 



QH 2 



(2) 



as the magnetic Reynolds number, and for small magnetic 
Prandtl number Pm = v/r}, the viscosity hardly influences 
the instability map. At a distance of 1 AU from a central 
mass with 1 Mq the disk thickness is about 3%, the density 
p of the gas 10 -10 g/cm 3 , and r\ ~ 4 • 10 15 cm 2 /s (because 
of the low temperature of the gas). Hence, the surprisingly 
high value of 0.1 gauss is obtained for the minimum B z . 

In order to excite nonaxisymmetric MRI modes one 
needs even higher magnetic fields. Kitchatinov & Riidiger 
(2010) showed in a linear theory that nonaxisymmetric 
modes only arise if Qa ~ 0.05 Q, with the Alfven frequency 
= Bz/^JJmTpH. For (2 = 2 • 10~ 7 at 1 AU this condition 
is fulfilled for B z > 0.3 gauss. Such strong fields at this dis- 
tance cannot be due to a stellar dipole field in the center. On 
the other hand, one should mention that the fossil magnetic 
fields found in meteorites are of just this order. One solution 
of this dilemma is to consider toroidal fields, which by the 
stretching action of differential rotation may often exceed 
the values of the poloidal field. For protoplanetary disks, 
however, the magnetic Reynolds number Rm only reaches 
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values of about 50, so that the resulting toroidal field does 
not dominate the original B z . 

The same is true for galaxies when the interstellar tur- 
bulence is driven by SN-explosions. The resulting magnetic 
Reynolds number is then also of order 10-100, so that the 
induced toroidal field does not dominate the poloidal one - 
as observed. 

A very different situation holds for the radiative cores 
of rotating stars. Due to the high conductivity of the hot 
plasma the magnetic Reynolds number is so large that even 
very slight differential rotation will produce toroidal fields 
which clearly dominate the fossil poloidal fields. It should 
then be sufficient to consider the stability of only the strong 
toroidal field under the presence of differential rotation. In 
the present paper we shall always apply a rotation law close 
to the profile of uniform specific angular momentum, i.e. 
Q oc R~ 2 . This rather steep profile is just beyond the regime 
of hydrodynamic centrifugal instability, also in spherical sys- 
tems (see Riidiger & Kitchatinov 1996). 

There are basic stability conditions in cylindrical ge- 
ometry. A flow is stable against axisymmetric perturbations 
if 



dR 



(R 2 n) 
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(Michael 1954). This condition turns into the well-known 
Rayleigh criterion for both B^ = and B^ oc R. The latter 
case means that an axial uniform electric current does not 
suppress Taylor vortices, which are due to too steep rotation 
profiles. Hence, rotation laws which are steeper than Q oc 
RT 2 will always excite axisymmetric rolls, even under the 
presence of very strong toroidal fields with B$ oc R. On the 
other hand, any unstable ^-profiles can be stabilized with 
a toroidal field of sufficiently large amplitude and suitable 
geometry. 

If nonaxisymmetric modes are considered then the nec- 
essary and sufficient condition for stability becomes 



dR 



{RB'l) < 



(4) 



(Tayler 1957, 1973, Vandakurov 1972), at least in the ab- 
sence of rotation. The fully general formulation including 
rotation is not known. Two questions immediately arise: i) 
is the field corresponding to a uniform current (_B^> oc R) 
really unstable, and ii) does the current-free field B^ ~ R" 1 
become unstable under the presence of (differential) rota- 
tion? Both questions can be answered. In the first case the 
resulting instability is called the Tayler instability (TI), and 
in the second case, where fields and rotation profiles that 
would each be stable on their own interact to yield insta- 
bility, the result is called the Azimuthal MagnetoRotational 
Instability (AMRI). Whereas the energy for the TI comes 
from the electric current, the energy source for AMRI is 
entirely from the differential rotation, and AMRI can only 
saturate at the expense of the rotational shear by producing 
high values of 'eddy' viscosity. This could therefore justify 
the high values of 'artificial' viscosity which Eggenberger, 
Montalban & Miglio (2012) introduced to explain the slow 
rotation of the stellar cores after the collapse towards their 
red giant stage. 

It is possible to demonstrate the mechanism which leads 
to AMRI with a simple dispersion relation for the stability 



of nonaxisymmetric perturbations in the form exp( r yt + kz + 
m(j)) with $t(-y) as the growth rate. It is enough to apply 
the case with Pm — s> and Q oc 1/R 2 (the Rayleigh limit) 
to the expressions given by Kirillov, Stefani & Fukumoto 
(2012) for nonaxisymmetric perturbations, which leads to 
the dispersion relation of second order 

Re* 2 ui 2 + Re* (03 + i&3)w + 04 + i&4 = 0. (5) 

Here we have used the quantity u) = j/VPm- Re with Re = 
Q/uk 2 as the growth rate in units of the angular velocity 
of rotation. The imaginary part of u describes an azimuthal 
drift of the nonaxisymmetric pattern (see below). 
The coefficients a and b are 



a s = 2(1 + (2 + m* 2 Ha*)), 63 = 2m* Re* 
04 = 1 + (4 + 2m* 2 ) Ha* 2 + m* 4 Ha* 4 — m* 2 



Re* 



&4 = 2Re m 
with the abbreviations 
Re* = a Re, Ha 



l + m* 2 Ha* 2 ) 



(6) 



— Ha, m* = — , (7) 
kR a 

where a = k z /k and Ha = B^/k ■ ^/Jlopvfj. 

The formal solution of ([5} for marginal instability leads 

to the two conditions 

64 



03(0304 + 0304) 



(8) 



Re* 03 ' 

The first of these relations gives the negative drift rate in 
units of the basic rotation, i.e. 



ojdr m*(l + m* Ha* 
~Q ~ ~ 1 + (2 + m* 2 )Ha* 



(9) 



The observable drift in the laboratory system, d<f>/dt = 
— cjdr/m, is thus positive (in the rotation direction) and with 
d(j)/dt ^ kQ jk z it is (for large Hartmann numbers) indepen- 
dent of the magnetic field, and also of the sign of the mode 
number m. 

The second condition in dsl leads to 



Re* 



3 Ha* 



(10) 



the latter relation holding for strong fields. For medium mag- 
netic field the relation changes to a linear one. Axisymmet- 
ric solutions do not exist. It is clear also that the mode with 
m = 1 is most easily excited: 

p 1 Ha2 (1U 
Re "2fcf^' 

valid again for large Ha. Solutions with large axial wave 
number - representing oblate cells - need much slower rota- 
tion to be excited. 

For the Elsasser number A = B^/ [lopnQ we obtain the 
relation 



A~ 2{k z R,f 



(12) 



where the influence of the microscopic viscosity completely 
disappears. 

With view on the experimental demonstration of AMRI 
we consider a Taylor-Couette set-up with two corotating 
cylinders of nearly perfectly conducting material. The cylin- 
ders confine an incompressible, conducting fluid under the 
presence of a toroidal magnetic field. It is clear that the ra- 
dial profiles of the rotating flow and the field between the 
cylinders are 
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Figure 1. The AMRI at the Rayleigh limit, fip = 0.25, for the two magnetic Prandtl numbers Pm = 10 -6 and Pm = 1. Left: the 
instability map, middle: the azimuthal drift of the solutions normalized with J7j n , right: the axial wave number multiplied with the gap 
width. 



U 4> = RQ = a sl R+^, B 4 , = a B R+^. (13) 

Let the ratio of the rotation rates of the cylinders be fin = 
^out/An and a similar expression for the field amplitudes, 
i.e. 



The ratio of the cylinder radii is n n = Ri n /R ut, where we 
will fix n n = 0.5. For such a container fin = 0.25 describes 
the Rayleigh limit for which R 2 Q ~ const ('shellular' ro- 
tation, an = in Eq. (|13[) L while fin — 0.5 describes the 
quasi-galactic rotation law U$ ~ const. Only these two ro- 
tation laws will be considered in detail here. 

For the magnetic field the expressions in ()13|) also con- 
tain two extrema. With &b = the toroidal field is due to 
a homogeneous current inside the fluid. For r- m = 0.5 the 
corresponding value is fis = 2. Such fields are unstable ac- 
cording to the criterion (Q against nonaxisymmetric pertur- 
bations. The existence of this kink-type 'Tayler' instability 
has recently been shown in the laboratory (Seilmayer et al. 
2012, Riidiger et al. 2012). 

The equations and boundary conditions are given by 
Riidiger et al. (2007) and will not be repeated here. We 
imagine the cylinders, which are unbounded in the axial 
direction, to be made of a perfectly conducting material. 
AMRI in the described container appears for fis = 0.5. The 
electric current only flows within the inner cylinder; no back- 
ground electric current flows through the fluid between the 
cylinders. 

For the rotation law with fin = 0.25 (the Rayleigh 
limit) the plots in Fig. [1] demonstrate the characteristic val- 
ues of the solutions. The critical Reynolds numbers (left), 
drift rates (middle) and wave numbers (right) are given for 
the magnetic Prandtl numbers 10 -6 and 1, which differ by 
many orders of magnitudes. The differences of the charac- 
teristic values given in Fig. [TJ however, remain very small. 
The AMRI at the Rayleigh limit can be more easily excited 
for Pm = 1 than for Pm — > 0, but the effect is weak com- 
pared with the variation of the magnetic Prandtl number. 
Obviously, the AMRI close to the Rayleigh limit scales for 
various Pm with the Reynolds number and with the Hart- 
mann number (Hollerbach et al. 2010). As an immediate 
consequence, the existence of the AMRI can be proven by 
laboratory experiments with liquid metals with their very 
low values of Pm -just as the helical MRI (Stefani et al. 
2009), but in strong contrast to the standard MRI which is 
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Figure 2. The normalized growth rate oj gr = 3(aj)/J?i n as func- 
tion of the toroidal background field, m = ±1. fin = 0.25, 
Re = 3000, Pm = 10" 6 . 



much harder to achieve (Riidiger & Hollerbach 2004). The 
results obtained with materials with low Pm are also repre- 
sentative of materials with Pm of order unity. Typical values 
for such possible experiments are Re ~ 1000 and Ha ~ 100. 

The meaning of the marginal instability curves is ex- 
plained by the growth rates, given in Fig. [2] for a fixed 
Reynolds number Re = 3000. The growth rate is positive 
between the two critical values, which in Fig. Q] (left) en- 
close a region of instability. The maximum growth rate of 
0.02 lies between the two limiting Hartmann numbers, and 
implies a growth time of about 7 rotation times. Compared 
with the standard MRI the AMRI is slower. 

As the left-hand plot of Fig. [T] also shows, the AMRI 
instability domain is always limited by two values of the 
Reynolds number or two values of the Hartmann number. 
For fixed rotation rate the magnetic field can be too weak 
or too strong, and conversely for fixed magnetic field the 
rotation can be too slow or too fast. 

The azimuthal drift rates (Fig. [I] middle) are always 
negative, hence the pattern of the instability wave drifts in 
the direction of the basic rotation (see above). A typical 
value of the drift (normalized with Oi n ) for small Pm is 0.25 
(marked in the plot), which hardly depends on the Hart- 
mann number. For both the Rayleigh limit and the rotation 
ratio fin = 0.25 the pattern essentially corotates with the 
outer cylinder. For stronger fields the drift is slightly slower. 
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Figure 3. The maximal normalized growth rate as a function 
of the Reynolds number for fin = 0.25. The Hartmann numbers 
which correspond to the maxima vary from about 200 to about 
500 from the left to the right axis of the plot. Pm = 10 -6 . 



The right-hand plot in Fig. [T] gives the axial wave num- 
bers normalized with the gap width d — R ou t — Rin- A wave 
number of n (marked) describes a cell with a circular geom- 
etry in the radial and the vertical coordinates. Precisely this 
cell form exists for the instability at the weak-field end of 
the instability domain. For stronger fields the cells become 
more oblate (not prolate!) but this is only true for the lower 
branches of the instability map. Along the upper branches 
the cells preserve their circular shape. 

One must ask how the growth rate behaves for increas- 
ingly rapid rotation. To answer this question one must deter- 
mine the maxima known from Fig.[2]for increasing Reynolds 
numbers. The results given in Fig. [3] clearly show the sat- 
uration of oj gr = Q(uj)/Qi n . The maximum growth rate for 
very rapid rotation is 0.045 J?i n , so that the shortest growth 
time of AMRI at the Rayleigh limit is 3.5 rotation times. 



2 QUASI-GALACTIC ROTATION 

For the quasi-galactic rotation law with \xn = 0.5, Fig. [4] 
shows the domain of instability in the S-Rm plane, for vari- 
ous magnetic Prandtl numbers. It is obvious that for Pm < 1 
the values scale with the magnetic Reynolds number Rm and 
the Lundquist number S, i.e. 



Rm = 



l?in d 



B in d 



(15) 



The unstable domain in the plane again has the characteris- 
tic form of tilted cones, so that minimal and maximal values 
of both S and Rm always exist for the onset of the instabil- 
ity. Hence, both the magnetic field and the rotation rate can 
be too weak or too strong for AMRI. One also finds that the 
instability curves for Pm — > converge in the S-Rm plane. 
For all Pm the ratio Oa/£2 lies between a low- field limit and 
a high-field limit, e.g., for Pm = 1, 



0.3<£<1. 



(16) 



It is thus clear that for AMRI the Alfven frequency of the 
toroidal magnetic field and the rotation rate must be of the 
same order. If the toroidal field is produced by differential 
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Figure 4. The instability map for the quasi-galactic rotation law 
fiQ = 0.5 in the S — Rm-plane. m = ±1. The curves are marked 
with their values of the magnetic Prandtl number. 



rotation acting on a poloidal fossil field, then the low-field 
limit in Eq. (|16[l plays the role of the onset condition for the 
axisymmetric instability. 

For given Rm and Pm the growth rate has been calcu- 
lated between the two limiting values S where it vanishes; it 
is maximal somewhere between the two limits. In Figs.[5]the 
ratio uj s 

Pm one finds a quasi-linear relation 



g r/J?m is plotted for various Rm and Pm. For fixed 
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or 



Rm 



(it 



with £ gr of order 10 



(17) 



(18) 



It varies from 1.5 ■ 10 for Pm = 1 



to 2.1 • 10~ 4 for Pm = 0.01, which suggests a very weak 
dependence of the growth rate cj gr on Pm. Hence, the growth 
time in units of the rotation time is r gr /r rot ~ 10 3 /Rm for 
small Rm. For smaller Rm the AMRI is rather slow, but the 
linear relation (|17p can only hold for small Rm. Figure 
shows that the growth rate slowly grows for smaller Pm but 
this effect is rather weak. 

The saturation of the normalized growth rate for Pm = 
1 is demonstrated by Fig. [6] For sufficiently rapid rotation 
the dependence of uj gI on the value of Rm disappears, so that 
always SJ(w)/i7j n ^ 0.14. The growth time, therefore, for the 
instability of the considered rotation law and for Pm = 1 can 
never be shorter than one rotation time. 



3 DESIGN OF AN EXPERIMENT 

With the given data it is easy to design an experiment for the 
realization of AMRI in the laboratory. For experiments on 
magnetically induced instabilities of the differential rotation 
beyond the Rayleigh limit it is natural to work with a Taylor- 
Couette container with a rotation law of fj,n = 0.26, which by 
itself would be hydrodynamically stable. For Re = 3000 Fig. 
[7] shows that the limit of marginal instability is reached for 
a Hartmann number of 85. Note that the instability region 
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Figure 5. The normalized growth rate tj gr = S(tAi)/J?i n as a function of the toroidal background field. The curves are marked with their 
values of the magnetic Reynolds number, m = ±1. fin = 0.5, Pm = 1 (left), Pm = 0.1 (middle), Pm = 0.01 (right). 
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Figure 6. The normalized growth rate aj gr = 3(tj)/!2i n for 
Hq = 0-5, optimized in the tilted instability cone for Pm = 1 
as a function of S (top) or Rm (bottom) . 



is slightly smaller than for the rotation law with fig — 0.25, 
which is understandable as the energy for AMRI is com- 
pletely provided by the energy of the shear. 

The relation 7 ax i s = 5-Rm-Bm connects the toroidal field 
amplitude Bi n at Ri n with the axial current inside the inner 
cylinder. /, R and B must be measured in ampere, cm and 
gauss. Hence, 



Ha= -- 



5 ^fiopurj 



(19) 



The radial size of the container does not occur in this re- 
lation between the Hartmann number and the current am- 
plitude. For the gallium alloy GalnSn the magnetic Prandtl 
number is 1.4 ■ 10~ 6 and the value of the square root in (|19p 
is 25.6. The resulting electric current for marginal instability 
is 10.9 kA. 

For a container with (say) Ri n = 4 cm and with the 
viscosity of GalnSn (y = 3.4 • 10~ 3 cm 2 /s) the Reynolds 
number Re = 3000 requires a rotation rate of the inner 
cylinder of only 0.1 Hz. Figure [5] shows that for the weakest 
fields the cell structure of the nonaxisymmetric pattern is 
circular and that the pattern nearly drifts with the same 
rotation rate as the outer cylinder. Both properties are very 
characteristic for the instability close to the marginal limit 
for weak magnetic fields - independent of the steepness [iq of 
the rotation law. Such an experiment is currently operating 
at the Helmholtz-Zentrum Dresden-Rossendorf. The results 
will be presented in a separate paper. 
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Figure 7. The normalized growth rate for an experiment with 
Hq = 0.26 as a function of the toroidal background field, m = ±1 
, Re = 3000, Pm = 10" 6 . There are no visible differences to the 
results for Pm = 0. 



0.0 F 
-0.1 : 

CT -0-2 : 

3 - -0.3 - 

2 V -0.4 - 

ot , , , , ] -0.5 6 



100 200 300 400 
HARTMANN NUMBER 



100 200 300 400 
HARTMANN NUMBER 



Figure 8. The normalized wave number (left) and drift rate 
(right) for fixed Reynolds number Re = 3000. The dotted lines 
give the wave number of n for a circular cell structure (left) and 
the drift rate of —0.25 for corotation of the pattern with the outer 
cylinder (right), m = ±1, jifi = 0.26, Pm = 10 -6 . 



4 KINETIC AND MAGNETIC ENERGY OF 
THE INSTABILITY 

We have seen that the modes for m = 1 and m = — 1, cor- 
responding to left and right spiraling modes (Hollerbach et 
al. 2010) are fully degenerate, i.e. they are excited at the 
same eigenvalues and possess the same wave numbers and 
drift rates. For the following nonlinear calculations the spec- 
tral cylindrical MHD code of Gellert, Riidiger & Fournier 
(2007) is used. The solution is expanded in azimuthal Fourier 
modes. The resulting meridional problems are solved with a 
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Figure 9. The nonlinear instability pattern of AMRI is a com- 
bination of the modes m = ±1 which drifts with a common drift 
rate. There are seven cells along the axial direction with the as- 
pect ratio r = 8. The cells are thus slightly elongated. Re = 500, 
Ha = 80, fia = 0.26, Pm = 0.1. 



Legendre spectral element method after Deville, Fischer & 
Mund (2002). The code, however, is only able to solve equa- 
tions with a minimum magnetic Prandtl number of 0.01. Fig. 
[9] shows the structure of the resulting wave, which consists 
of both m = ±1 components, and drifts together with the 
outer cylinder. The almost circular meridional cell structure 
also confirms the results of the linear theory. There are seven 
cells along the vertical axis with its eight length units. 

Besides the geometric structure, the nonlinear code also 
provides the behavior of the equilibrated energies. The ques- 
tions are how strong the magnetic energy is (i) in units of the 
kinetic energy, and (ii) in units of the energy of the toroidal 
background field. The maximal amplitude of the radial field 
component in units of Bi n in Fig. [9] is rather small. This 
fact, however, is only due to the rather slow rotation of this 
example. 

Figure [10] displays the energies of several instability 
models for magnetic Prandtl numbers between 0.05 and 1.0. 
Both the magnetic and the kinetic energy are integrated 
and normalized with the energy of the toroidal background 
field. All energies are integrals over the entire container. The 
first finding concerns the quantity q = {b 2 )/{Bf n }, which 
for driven turbulence in the high-conductivity limit simply 
scales with Rm (Brauer & Krause 1974). The same is true 
for this quantity in the case of AMRI. From Fig. QJj] the 
relation is 



(Rm — Rm c 



(20) 



with e mag — 6 • 10 -4 . It should thus be possible that the 
energy of the magnetic fluctuations reaches the order of the 
magnetic background field. The relation ()20[) can also be 
written as q oc £ ma gI2 A»V The saturation of the relation (|20p 
for large magnetic Reynolds numbers, however, is not yet 
known. A possible form of the relation might be q oc Q\j Q 1 . 
Figure [10] also shows that the magnetic energy always 
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Figure 10. The magnetic (solid) and the kinetic (dashed) ener- 
gies of AMRI for various magnetic Prandtl numbers vs. the mag- 
netic Reynolds numbers. Re = 1000, Ha = 75 . . . 200, ixq = 0.26, 
Pm = 0.05. . . 1. 
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Figure 11. The kinetic energy (solid) and its part in the radial 
direction (dashed) of the AMRI perturbations for various Pm vs. 
the modified magnetic Reynolds number Rm. Ha = 75 . . . 200, 
HQ = 0.26. 



dominates the kinetic energy of the fluctuations. One may 
speculate that the factor of four which arises here also ap- 
pears in the ratio of the resulting eddy viscosity and the 
chemical diflusivity (see below). According to Fig. 1111 the 
kinetic energy for various magnetic Prandtl numbers vs. the 
modified Reynolds number Rm = Rm/vPm = Qi n Rf n /f] 
with fj — yfvr\ provides the linear relation 



B? 

in 

with £kin 



Skin (Rm 



1.4- 10" 



- Rm C rit) 
Hence, 



Mrms -0.012 A 5 



with the modified Elsasser number 

A- Bl 



(21) 



(22) 



(23) 



which does not depend on the size of the container. The 
appearance of the molecular viscosity in the expression of 
the kinetic energy might have strong consequences for the 
value of the kinetic energy. As the ratio of the magnetic 
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Figure 12. The proxy (125 1 of the turbulent Schmidt number for 
various magnetic Prandtl numbers and for fiQ = 0.26 does not 
depend on the magnetic Reynolds number but decreases as \/Pm. 



and kinetic energies scales with VPm it should be possible, 
for small magnetic Prandtl numbers, that the kinetic energy 
even dominates the magnetic energy. For not too small Pm, 
however, it is clear that the AMRI provides higher values of 
the viscosity than for the diffusivity of chemicals. The an- 
gular momentum transport is dominated by the magnetic 
energy while magnetic fluctuations do not contribute to the 
diffusion of thermal energy and/or the mixing of chemicals 
(Vainshtein & Kichatinov 1983). The different scaling of the 
magnetic and kinetic energies has consequences for their ra- 
tio e = (b 2 ) I ^op(u 2 ) , which obviously scales as VPm. The 
magnetic energy only dominates for magnetic Prandtl num- 
bers of order unity or larger. Consequently, the turbulent 
Schmidt number 
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should also scale with VPm. Here en is formed with the 
radial velocity components (also given in Fig. (jlll) . i.e. 



£r = 



{b 2 } 



(25) 



as only the radial turbulent intensity is responsible for the 
radial mixing of chemicals. Figure [12] demonstrates that in- 
deed the turbulent Schmidt number on the basis of AMRI 
strongly exceeds unity, but decreases with VPm for small 
Pm. Note, however, that for very small magnetic Prandtl 
number it cannot become smaller than 0.4 (Yousef et al. 
2003). In media with small magnetic Prandtl number the 
angular momentum transport and the mixing of chemicals 
by AMRI are of the same intensity. Both effects only stop 
when the regime of rigid rotation is reached. As long as the 
differential rotation is strong enough the mixing of chemi- 
cals is also strong. It might be, however, that the inclusion 
of the 'negative' buoyancy by the stable stratification of real 
stellar cores only reduces the mixing rather than the angu- 
lar momentum transport and the Schmidt number is not 
reduced too much by the small magnetic Prandtl number 
(Kitchatinov & Brandenburg 2012). 



